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ABSTRACT 

The  angle  ♦  which  the  free  boundary  of  an  extreme  wave  makes  with  the 
horizontal  is  the  solution  of  a  singular,  non-linear  integral  equation  that 
does  not  fit  (as  far  as  we  know)  into  the  theory  of  compact  operators  on 
Banach  spaces.  It  has  been  proved  only  recently  that  solutions  exist  and  that 
(as  Stokes  suggested  in  1880)  these  solutions  represent  waves  with  sharp 
crests  of  included  angle  2*/ 3.  In  this  paper  we  use  the  integral  equation, 
known  properties  of  solutions  and  the  technique  of  the  Mellin  transform  to 

obtain  the  asymptotic  expansion 

_  Jj  V_  Pk 

♦(b)  ■  r  +  V  asn  +  o(s  *)  as  s  +  0  ,  (*) 

t>  _  n 

n*1 

to  arbitrary  order;  the  co-ordinate  8  is  related  to  distance  from  the  crest 
as  measured  by  the  velocity  potential  rather  than  by  length.  The  first  few 
(and  probably  all)  of  the  exponents  Mn  are  transcendental  numbers.  We  are 
unable  to  evaluate  the  coefficients  an  explicitly,  but  define  some  in  terms 
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of  global  properties  of  4*  and  the  others  in  terms  of  earlier 
coefficients.  The  derivation  of  (*)  includes  an  assumption  about  a  question 


in  number  theory;  if  that  assumption  should  be  false,  logarithmic  terms 


P  ra 

n. .  .  n § j 

a  .s  (log  s )  , 

n*  j 


with  in  .  a  positive  integer,  would  enter  the  series  at 
n*D 


very  large  values  of  n. 

Our  results  confirm  the  heuristic  calculations  of  Grant  [4]  and  Norman 
[9] ,  and,  in  effect,  remove  a  tentative  element  from  certain  statements  in 
those  papers. 


AMS  (MOS)  Subject  Classifications:  76B15,  45G05. 

Key  Words:  water  waves,  non-linear  integral  equations,  asymptotic  analysis. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  concerns  waves  of  permanent  form  on  the  free  surface  of  an 
ideal  liquid  which  is  in  two-dimensional,  irrotational  motion  under  the  action 
of  gravity.  We  consider  only  extreme  waves,  often  called  'waves  of  greatest 
height' t  each  of  these  is  the  end-member  of  a  one-parameter  family  of  waves, 
and  is  distinguished  from  other  'smaller'  members  of  the  family  by  a  sharp 
crest «  aa  in  Figure- -ila^  Although  this  corner  is  physically  unrealistic, 
oceanographers  have  given  such  idealized,  extreme  waves  a  great  deal  of 
attention  since  Stokes  postulated  their  existence  in  1880.  (One  reason  may  be 
the  physical  importance  of  the  smaller  waves,  and  that  scientists  like  to 
interpolate. ) 

The  present  paper  is  a  contribution  to  the  strict  mathematical  theory  of 
extreme  waves,  which  has  emerged  only  since  1978.  We  derive  rigorously  an 
asymptotic  series  that  describes  the  flow  near  the  crest.  T^iis  confirms  and 
sharpens  certain  earlier  exploratory  results  due  to  Grant  -H^>and  Norman .  . 

The  series  should  play  a  useful  part  in  numerical  computation  of  extreme 
waves. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  authors  of  this  report. 
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OH  THE  BEHAVIOR  HEAR  THE  CREST  OF  HAVES  OF  EXTREME  FORM 


C.  J.  Amick1  And  t.  E.  Fraenkel2 


1.  Introduction 

This  pspsr  concerns  gravity  waves ,  of  permanent  and  extreme  form,  on  the  free  surface 
of  an  ideal  liquid,  the  floe  being  two-dimensional,  irrotational  and  in  a  vertical  plane. 

By  a  wave  of  extreme  form  we  mean  one  that  is  the  'largest'  member  of  a  one-parameter 
family  of  such  waves,  and  is  characterised  by  a  sharp  crest  of  included  angle  2x/3,  as 
shown  in  Figure  1(a).  The  existence  of  such  waves  was  conjectured  by  Stokes  in  1880  and 
has  been  proved  recently.  (For  a  fuller  account,  see  the  Introduction  to  [1]»  for  back¬ 
ground  material,  see  [8].) 

letting  4  denote  the  local  wave  angle  (that  is,  the  angle  which  the  free  boundary 
makes  with  the  horisontalf  *(s)  •  tan  V' (x)  in  the  notation  of  Figure  1),  we  shall  seek 
the  behavior  as  s  ♦  0  of  a  solution  4  of  the  equation 


♦o-jj;  K(..t)  y<|>  nt>  dt,  o<.<*  , 

J0  v  ♦ 


(1.1) 


where 


and 


tan  2  *  +  tan  2  t 


!C(S,t)  -  ^  log- - , - r~ - 

I tan  -  s  -  tan  -  t| 


v(t)  -  -j  (cos2  |  t  +  b  sin2  -j  t) 


for  son*  b  8  [0,1]  which  we  regard  as  fixed  henceforth.  Here  jg  v  sin  ♦  stands  for 
v(u)sin  S(u)duf  such  abbreviations  will  be  used  throughout  the  paper.  Equation  (1.1) 
refers  only  to  waves  of  extreme  form  the  end-value  of  a  parameter  has  already  been  chosen. 
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The  significance  of  the  weighting  function  v  is  that  v(s)  «  -C  34/3s,  where  C  is  a 
positive  constant  and  4>  denotes  the  velocity  potential;  the  constant  b  distinguishes 
the  families  mentioned  earlier:  b  »  0  corresponds  to  solitary  waves  (of  infinite  wave¬ 
length  and  in  liquid  of  finite  depth),  b  e  (0, 1)  to  periodic  waves  in  liquid  of  finite 
depth,  and  b  »  1  to  periodic  waves  in  liquid  of  infinite  depth. 

By  a  solution  of  (1.1)  we  mean  a  function  t  satisfying  the  equation  pointwise  and 


such  that 


0  <  ♦(s)  <  —  on  (0,»),  $(s)  -  0,  ♦  is  real-analytic 


on  ((>,»],  and  $(e>  ♦  —  as  a  *  0 

o 


Such  solutions  are  now  known  to  exist  ((1],  [2],  [3],  [7],  [13]). 

There  have  been  many  proposals  for  calculating  the  shape  of  the  free  boundary  of  an 
extreme  wave  by  a  combination  of  analytical  and  numerical  methods  ( long  lists  of  references 
are  given  in  [5]  and  [14]).  Such  calculations  were  given  a  new  direction  by  Grant  [4],  who 
sought  the  second  term  of  an  asysiptotic  expansion  of  z(x)  for  x  *  0  (here  z  «  x  +  iy 
denotes  the  complex  co-ordinate  in  the  physical  plane,  and  x  *  *  +  l*  the  complex 
potential),  the  first  term  being  given  by  Stokes's  comer  flow.  Grant  pointed  out  that  the 
second  term  must  have  an  exponent  that  is  irrational  and  ’probably  transcendental*;  he 
concluded  that  ’the  structure  near  the  comer  is  considerably  more  complicated  than  has 
been  assumed  in  the  past*.  Norman  [9]  contemplated  terms  beyond  the  two  considered  by 
Grant;  inferred  the  nature  of  all  the  exponents;  introduced  the  assumption  that  the  numbers 
8^,  defined  after  (1.3)  below,  are  linearly  independent  over  the  rationals;  and 
established  certain  relationships  between  the  coefficients  of  the  series.  (However,  it 
seems  that  no  coefficient  after  the  first  can  be  calculated  by  a  merely  local  analysis.) 
Longuet-Higgins  and  Pox  [5],  [6]  used  the  exponents  and  functions  arising  in  the  work  of 
Grant  and  Norman  as  one  part  of  their  theory  of  waves  near  an  extreme  we  (the  crest  being 
smooth  but  of  very  large  curvature ) . 

All  these  calculations  were  heuristic  rather  than  rigorous;  in  particular,  Gkant  [4] 


and  Norman  [9]  were  concerned  only  to  find  analytic  functions  s(X>  that  satisfy  the  non- 
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linear  boundary  condition  of  the  problem  in  an  asymptotic  sense  as  ♦  ♦  0.  In  this  paper 


we  proceed  from  (1.1)  and  (1.2)  to  prove  that,  subject  to  an  assumption  about  a  question  in 
number  theory  to  be  explained  presently, 

"  u 

$(s)  ~r+  )  as"  as  s+0  .  (1.3) 

6  .  n 

n»1 

The  exponents  u  depend  on  the  positive  zeros  B  ,  B  ,  B  ,...,  with  B  <  8  ,  of 

n  i  2  j  j  i 

/3(1+B>  -  tan(Bx/2)»  the  6^  are  essentially  the  numbers  discovered  by  Grant  [4],  and 
6 1  ~  0.8027,  B2  ~  2.9066.  As  Grant  thought  probable,  and  as  we  prove  in  the  Appendix,  each 
B  is  a  transcendental  number.  Each  exponent  U  is  a  finite  linear  combination,  with 

j  n 

positive  integer  coefficients,  of  numbers  in  the  set  {2,6 ^ ,B2» . . .} ,  and  contains  at  least 

on#  B^ .  The  u  are  ordered  by  u  <  v  (as  the  symbol  ~  of  asymptoticity  in  (1.3) 
j  n  n  n+ 1 

implies);  the  first  few  are  B1»  28^,  36^  8^  +  2,  B2,  46^,  26^  +  2,  B?  +  B2»  56^, 

38 1  +  2,  26 1  +  B2,...  . 

The  assumption  made  in  the  derivation  of  (1.3)  is  stated  precisely  in  Section  4;  here 

we  remark  that  it  certainly  holds  if  the  set  { 1,$1 ,B2, . . .}  is  linearly  independent  over 

the  rationals.  Moreover,  numerical  calculation  indicates  that  the  assumption  is  true  for 

the  first  hundred  of  the  slightly  larger  set  of  exponents  arising  in  the  derivation  of 

(1.3).  If  the  assumption  should  be  false,  then  our  method  would  still  be  applicable,  but 

U  m 

logarithmic  terms  a  .a  "(log  s)  n'3,  with  m_  ^  a  positive  integer,  would  enter  the 
n  #  j  n  *  j 

series  at  large  values  of  n. 

He  cannot  evaluate  the  coefficients  a_  in  (1.3)  when  ti  e  (B  ,,B_,B  , . . . } ;  we 

n  n  i  i  s 

define  such  coefficients  by  integrals  involving  the  global  behavior  of  ♦ .  When 
U  4  {Bj»Bj#Bj»...},  the  corresponding  coefficisnt  an  is  determined  by  the  previous 
coefficients  •i'***'an-ii  this  agrees  with  Norman's  results  in  [9). 


The  expansion  (1.3)  can  be  transformed  and  integrated  to  yield 

•  l> 

z(X) - ig”1/3(|  1X)2/3(1  +  ).  bn ( ix )  ")  as  x  ♦  0  » 

n»  1 


(1.4) 


where  the  constant  g  is  the  gravitational  acceleration,  arg(ix)  e  I-*/2,*/2] ,  and  the 

coefficients  bn  are  real.  If  the  set  f 1,61 ,f*2, • . .)  is  linearly  independent  over  the 

rationals,  then  b_  «  0  whenever  the  linear  combination  defining  u  contains  a  multiple 
*»  n 
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of  2.  The  expansion  (1.4)  is  than  of  the  form  proposed  by  Norman  (9),  if  we  interpret 
liberally  certain  tentative  remarks  in  that  paper  (for  example,  that  it  is  'possible  to 
consider  solutions  ...  corresponding  to  combinations  of  terms  from  several  roots'  of 
✓  3(1+8)  «  tan(  Bw/2) ) . 

The  plan  of  the  paper  is  as  follows.  We  begin  Section  2  by  making  the  transformation 

C  «  tan  ~  s,  ♦  (?)  “  4(2  tan  *5)  in  order  to  obtain  a  kernel  k(C,n>  that  is  simpler  than 

K(s,t).  (In  effect,  we  map  the  unit  disk  in  the  plane  of  C  “  peiS  onto  the  half-plane 

{(i>  ■  a  +  iC  s  c  <  0}  by  the  conformal  transformation  a>  -  (C-1)/(C  +  1)  • )  Since  4(8)  • 

♦(tan  s),  an  asymptotic  expansion  of  1<(C)  for  C  +  0  yields  one  of  $(s)  for 

s  +  0.  The  next  step  is  less  obvious;  we  cast  the  integral  equation  (2.1)  for  ♦(£)  into 

a  form,  (2.10),  that  contains  an  elaborate  non-linearity  but  has  the  virtue  of  allowing  us 

to  construct  the  expansion  of  ♦(£)  by  an  inductive  process.  To  begin  this  process,  we 

a  13 

show  that  ♦(?)  -  */6  -  0(C  )  for  some  a  e  (-,  -) . 

4  4 

In  Section  3  we  combine  this  preliminary  estimate  with  the  use  of  the  Mellin  transform 
to  show  first  that 

0 

♦(C)  -  |  +  AtC  1  +  o(Cl>  as  5  +  o  (1.5) 

for  any  t  e  (8^1),  and  then  that  this  result  implies  the  improved  approximation 

8  28 

♦  (C)  -  l  +  A,C  1  +  A25  1  +  o(C">  as  C  *  0  0.6) 

for  any  m  e  (28^,8^  +  l)i  here  Aj  is  a  known  function  of  Aj. 

The  step  from  (1.5)  to  (1.6)  points  the  way  to  the  long  inductive  proof,  in  Section  4, 
of  the  main  results  of  the  paper,  which  appear  in  Theorem  4.5  and  Corollaries  4.6  and 
4.7.  The  Appendix  concerns  some  properties  of  the  numbers  8^. 
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2.  Preliminary  transformations  and  estimates 

Under  the  transformation  £  »  tan  s,  n  *  tan  -j  t  and  ♦  (£)  *  $(2  tan  \ )  =  <t>  (  s ) , 
equation  (1.1)  becomes 

♦  <C)  -4  j!  *(C.n)  dn.  o<£<-.  (2.1) 

30  Jp  w  sin  4> 

where 

k(£,n)  ■  ^  log  -j and  w(n)  -  ( 1  ■*•  n2)  ^2  ( 1  +  bn2)  ^2  • 

By  a  solution  of  (2.1)  we  mean  a  function  ♦  satisfying  (2.1)  pointwise  and  such  that 

0  <  ♦  (£)  <  ~  on  (0,-),  ♦  (£)  ■  0(£  1 )  as  £♦«•,♦  is  real-analytic  j 

)  (2.2) 

on  (0,»),  and  ♦  (£)•»•£  as  £  ■►  0  .  J 

It  is  to  be  understood  henceforth  that  £  e  ( 0  ,•) .  Occasionally  we  set  £  “  0 ,  with  the 
implication  that  ♦  (<))  «  */6  and  (i  e  C[0,«)  (even  though  in  the  original  problem  t|>  is 
an  odd  function  on  R\{0}). 

Combining  (2.1)  and  the  formula  ([1],  p.  197) 


T  r0  MC,n)  ^dn  -  J  , 


we  obtain 


♦(C)  “  g  *  "J  J*  MC.n)  ~  log{^  Jg  w  sin  ♦)dn  , 


for  any  constant  C  >  0.  Define 


7(C)  -  ♦(€)  -  (ET)(n)  -  2w(n)  sin{|-  +  y(n)}  ; 


7(C)  -  |  Jg  M£,n)  ~  log{^  jj]  EYldn  .  (2.4) 

Our  next  transformation  of  the  equation  is  more  elaborate;  it  involves  the  integrated 
kernel  (used  extensively  in  [1]) 

2  2 

q(C,h)  -  Jg  k(t,n)dt  “  ~  {c  log  |^|  +  n  log  ^  ~n  , 


qn(£,TD  -  7  log  i5-“nJ 

n 

and  the  non-linear  operator  F  defined  by 


°  "  ■  • 


. *1  *4  •  *  *  .»■  _•>  .*>  »  %  •*  *  * 


(2.5a) 


<FY)(n> 


log{n  C  EYl 


+ 


n(Et)  (n) 


i 


-  ~  In  log{^  /q  ®y)]  •  (2.5b) 

The  ‘unction  F*Y  is  important  throughout  the  paper  j  we  note  its  behavior  for  n  ♦  0  and 
for  n  ♦  •.  By  the  definition  (2.3), 

(FY)(n)  -  i  +  /i  Y(n)  +  o<Y(n)2  +  n2)  as  n  +  o  ,  (2.6) 

and  it  follows  froa  either  fora  of  (2.5)  that 

<rr)(n)  -  /3  Y(n)  +  o(y(i)2  +  n2)  as  n  ♦  o  .  (2.7) 


He  claia  that 


,  1  ,n 

o  <  /  by  <  J0  rr  <  const.,  i  <  n  <  -  , 
where  the  constant  depends  on  Y*  The  lower  bound  follows  froa  (2.6),  in 
as  n  ♦  0,  and  because  (EY)(?)  >  0  on  (0,»)  by  (2.2).  For  the  upper 

ij  EY  <  2  Jj  *(h)(1+h2)  /zdn  -  /’  4(t>  sec  “  dt  <  • 

by  [3],  p>  657,  if  be  (0,1],  and  by  [2],  Theorem  4.7(a),  for  the  case 

solitary  wave.  Also,  n(BY)(n)  is  bounded  because  w(n)  -  0(n  ’)  as  0 

(2.5a)  and  (2.8)  show  that 


(2.8) 

which  Y(h)  ♦  0 
bound  we  have 

b  »  0  of  a 
♦  <*  i  hence 


!(FY)(n)|  <  const,  log  n,  2  <  n  <  -  . 


(2.9) 


LEMMA  2.1.  If  W/6  ♦  Y  is  a  solution  of  (2.1),  then 

Y(C)  -  j  ip  r<e,n)(FY)(n)dh,  0<C<-  , 


(2.10) 


where 


r(C,n)  -  -  £  qn(C,n> 


Proof,  with  the  notation  (2.3),  equation  (2.1)  becoaes 


|  ♦  T«)  -  i  io  k<5'n>  dn  • 

Jo  W 


(2.11) 


Differentiation  of  this  with  respect  to  5  is  legitiaate  if  the  resulting  Integral  is 
written  as  a  Cauchy  principal  value;  noting  that  ?h^(C,h)  “  ,h) ,  and  integrating  by 

parts,  we  obtain 


CY’(C) 


it i  i“i* 

Jo  ** 


(2.12) 
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where  the  integral  exists  by  the  properties  {2.2)  of  \|i  =  n/6  +  Y  and  because 


«  0(n  )  as  n  *  0  (see  (7)).  We  add  (2.4)  and  (2.12),  refer  to  the  definition 
(2.5a)  of  FT,  and  integrate  with  respect  to  £;  there  results 

€y(5)  =  |  /”  q(C,h)  ( FY ) ( h ) dn  . 

We  may  integrate  by  parts  because  (for  fixed  £  e  (0,»),  as  elsewhere)  q(£,n)  is 
0(n  log  1/n)  as  n  *  0,  and  is  0(n  )  as  n  ♦  and  because  we  have  the  estimates 

(2.7)  and  (2.9)  for  FY.  Accordingly, 

£y(£)  =•  -  -j  /“  9n(E.h)(FY)(n)dn  ' 

and  this  is  (2. 10). 

LEMMA  2.2.  If_  ir/6  +  Y  is  a  solution  of  (2.1),  then  there  exists  an  exponent 
a  @  (■“,■”)  such  that  Y(?)  “  0(5°)  as  £  *  0. 

2  4  ..  —  — 

Proof .  in  this  proof  we  abbreviate  (FYHh)  to  F^(n),  and  similarly  for  other 
functions  of  the  same  kind. 

(i)  For  any  c  e  (0,1],  let  £  e  <0,c]  and  rewrite  (2.10)  as 

Y<£)  -  |  Jg  r(£,n)F^(n)dn  +  R  (£,c),  0  <  £  <  c  ,  (2.13a) 

where,  in  view  of  (2.7)  and  (2.9), 

I Ry< ^ ,c) I  “  1^  /"  r(C,n)FY(n)dn| 

2  2 

<  const. {/^  r|log  ^  t n  ^  I dT1  +  /,  ^2  l0g  n 

c  n  2  n 

<  const. {/,£  log  ~ —  du  +  £  /“  l0g2n  dn> 

u -i  n 

<  const.  E/c  ,  (2.13b) 

where  the  constant  depends  on  Y  but  is  independent  of  c.  It  is  natural  to  define 

G  (n)  -  2  sin{£  +  Y(n)}  -  1  » 

then  G^(h)  ~  /3  Y(h)  as  n  +  0  and  Y(n)  +  0,  and,  by  the  definition  (2.3), 

E^(n)  «  w(n)d  +  <3^(0)}  »  1  +  G^(n>  +  0(n2>  as  n  +  o  . 

Moreover,  we  can  so  define  U  ,  V  and  w  that 


n 


i  n  i 

L  g  ♦  om  ,  , 

J  0  y 


{i  %r0cY  +  °<n  »"  ■ 1 

Jn  Ev 


n 

'!  s 


v  (n) 


1o9<^  /!  V  “  log{l  !oGy  +  0(n2,)  =  “V"  !\ 3  gy  +  0<n2'  ' 

g  (n)  »  cos  r(n)  -  i  +  /I  sin  Y(n)  ■/In  (n)Y(n)  ; 


then 


G^n)  *  0,  u^ln)  ♦  1,  v^cn)  *  i  and  w^(n)  ♦  i  as  n  ♦  o  , 


and 


rY(n)  -  (v^(n)  -  u  (n)  -  GY(n>uY(n)}  ^|Jg^  GY<n)  +  «y(n ) 

-  /3  {v^(n>  -  u^(n)  -  G^(n)u^{n>>  ^  /nQ  w^r  +  /3  w^nmn)  +  6^(n) 


2 

where  fi^(n)  ■  o(n  )  as  n  ♦  0. 


Finally,  define  a  linear  operator  L  on  the  space  C[0,c)  by 
( Lf )  ( £ )  -  j:  ]CQ  r(£,n)t{vY(n)-UY(n>-GT(n)UY(n)}  £  w^f  +  WY(n)f(n)]dn 
and  the  integral  equation  (2.13)  becomes 


Y  (£)  -  (LY) (E>  +  S^(£,c),  0  <  E  <  c 


where 


because 


st<E,c)  »  P^tE.c)  +  J  /q  r(£,n)8^(n)dn  »  o(£/c)  , 

2 

|/°  r(£,h)fi  (n)dn|  <  const.  £2  f^'\loq  ^1:"|u2du  < 
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const .  cE 


(li)  For  every  e  8  (0,r],  choose  c  ■  c(e)  so  small  that 

4 


| v^(n)  -  u  (n)  -  GY(n)UY(n)j  +  |i  -  »Y(n)|  <  e  for  all  n  e  [0,c(e>] 


this  is  possible  by  (2.14).  For  any  a  8  (0,1)  and  any  e  e  {0,-],  define  the 


apace 


X  -  X  (c( € ) )  -  {f  e  C[0,c(e)]  S  Ifl  <  •}  t 

ft  tt  a 


where 


,f,«  “  sup0<£<c(e)  • 

Consider  the  linear  operator  L  defined  by  (2.15);  we  wish  to  show  that  L 


into  itself,  for  c  sufficiently  small  and  for  some  a  >  -,  and  that 


(2. 14; 


(2.15) 

(2.16a; 

(2.16b) 

(2.17) 

I 

Banach 


maps 


as  follows. 


ILI ,  5  sup  ,  ,  I  Lf  I  /«fl  <  1  . 

<a)  Hfex\{0>  a  a 
a 

It  is  easy  to  show  that  Lf  e  C[0,c(e)]  when  f  e  X  ,  and  we  estimate  *LfA 

ci  a 

|(Lt)<e)|  <  ~  J°(e>  |r(C,ni|  /q  I  f  I  +  (i+€»|f|ldn 

<  M+e)ifla  ig(e)|r<5#n)  |(-“  n°  +  na>dn 


<  1«H1  +  £->««  j“  7  | log  -Ll%1^-|n“ 


/3w 


l+o  a  J  0  C 


dn 


~~T  n+e)<1  +  5°  /"  (log  | u°du  , 


/3» 


where  0  <  £  <  c(e).  It  follows  that 


where 


lL,(a)  <  (1+e)<1  +  — )u<a)  , 


M(a)  -  —7-  /“  | log  l|uadu 

/3« 


1  ,/«•  .  u2  o.  .  .  1I//2.  1-u2  a.  , 

“  ”  Un  log  - u  du  +  2  /  log  — =-  u  du}  . 

/3n  u  | 1-u* | 

The  flret  of  these  two  Integrals  is  evaluated  In  Section  3,  see  ( 3. 1  let ) »  In  the  second,  we 


set  u  ■  1/x  and  Integrate  by  parts >  then 


1  ,  v  os  4  ,•  dx  , 


(2.18) 


To  obtain  a  simple  majorant,  replace  x  by  1  In  the  last  Integral,  which  is  then  easily 
evaluated)  thus 


U(a)  <  - - -  {tan  ~  ~  l°g(1  +  /2)}  «  T(a>,  say. 

✓3(1+0)  2  * 

Now  !(■“)  <  0.82,  t(4)  >  1.16  and  T'(o)  >  0  on  (0,1))  hence  there  exists  a  number 

2  4 

o  6  (r#4)  such  that  y(o)  <  T(a)  <  r(a)  <  1  for  all  a  e  [0,aj.  If  we  choose  e 

2  4 

sufficiently  small,  then  ILI^  <  1  for  all  a  e  [0,a]. 

(iii)  With  e  and  c  «  c(e)  now  fixed,  the  estimate  (2.16b)  of  S I  may  be  written 

S  (C)  •  0(C)  as  C  ♦  0.  Hence  S„  e  X  for  all  a  8  (0,1),  and  so  the  equation 
7  7  a 

f (C)  «  ( Lf ) ( C )  +  S  (C)  has  a  unique  solution  f  e  X  for  each  a  e  fO,a]  .  Since 
7  a  a 
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f  for  all  a  e  [0,a] 


X  c  X  ,  these  solutions  are  Identical;  in  other  words,  f 
a  0  <» 

because  each  f  e  X„.  We  know  that  Y  e  X  and  satisfies  (2.16a),  whence  Y  ■  f.  @  X 

(J  0  0  U  U 

for  all  a  e  [0,a].  in  particular,  Y  e  X_,  so  that  a  may  be  chosen  as  the  exponent  in 

a 

the  statement  of  the  lemma. 

Remark .  We  shall  see  from  Lemmas  3.1  and  3.2  that  the  best  possible  exponent  is  the 
number  6^  introduced  after  (1.3),  so  that 

S  x 

1  1 
-  tan  — r—  «  1  . 

1+fl,) 

(Bt  is  best  possible  if  A,  ^  0  in  Lemma  3.2.)  This  corresponds  to  only  the  first  term 
of  y(a)  in  (2.18),  so  that  precise  treatment  of  the  second  term  is  not  worth  while. 


-10- 


*.  r. 

•  V*  »  * 

*  %  • 

•  V  WVv  V 


3.  The  asymptotic  expansion  of  i|i(5)  to  three  terns 

Our  observation  in  (2.7)  that  <FY)(n)  ~  /3  7(h)  as  n  ♦  0  suggests  that  the  integral 
equation  (2.10)  be  written 


where 


Y(f.) 


r(5,n)Y(n>dn  +  p^(5). 


o  <  5  <  i  , 


(3.1) 


v*>  ■  3  /!  r(5.n){(FY)<n)  -  a  7(n)}dn  +  |  j;  r<5,n>(FY>(n>dn  .  (3.2) 

From  (2.7)  and  Lemma  2.2  we  have 

<FY)(n)  -  /3  y(n)  -  o(n2a>  as  n  ♦  0,  with  a  e  <■■;»”>  , 

2  4 

and  it  follows  (by  an  estimate  like  (2.17)  for  the  first  integral  in  (3.2),  and  an  estimate 
like  (2.13b)  for  the  second)  that  p^(5)  ”  0(5)  as  5  +  0. 

All  the  results  to  come  will  be  consequences  of  (3.1)  and  (3.2).  Our  plan  is  to  boot- 

(X 

strap  from  some  expansion  of  <i(5)  for  5  +  0  (at  present,  n/6  +  0(5  ))  to  an  estimate 
of  Py(5)  (at  present,  0(5)),  and  then  to  derive  a  more  complete  expansion  of  1M5 ) 
from  (3.1),  regarded  as  a  linear  aquation  for  Y  in  which  p^  is  'known'.  The  proof  by 
induction  in  8ectlon  4  will  follow  this  program  to  exhaustion)  without  the  more  explicit 
first  steps  in  this  section,  it  would  probably  be  incomprehensible. 

The  Mel tin  transform  ([12],  p.  7)  will  be  our  main  tool  for  analysis  of  (3.1).  Let 
f  b  piecewise  continuous  on  ( 0 ,“) ,  define  4(x,s)  -  xB_1f(x)  for  x  >  0,  where 
s  -  o+it  e  C,  and  assusw  that  4(*,s)  €  L^O,*)  for  a  <  o  <  b.  Then,  for  s  in  this 
strip,  we  define  the  Mellin  transform  f  of  f  by 

f(s)  -  Jp  x"_1f(x)dx  i  (3.3) 

we  also  write  this  as  f(s)  ™  f(x).  Note  that  f  is  analytic  in  the  strip  Be  s  e  (a,b). 
The  Inversion  formula  ([12],  p.  46)  is 

f(x)  "  *  "  *(■)<*■,  a  <  c  <  b  ,  (3.4) 

at  points  of  continuity  of  f.  The  product  formula  ([12),  p.  54)  is 

f (-s+A)g(s+B) —  x  J0  u  f(u)g(xu)du  ,  (3.5) 


provided  that  f(-s+A)  and  g(s+B)  have  a  common  strip  of  convergence.  Finally,  we 
record  the  following  property  ([121,  P*  118)  of  the  Mellin  transform. 


-12' 


LEMMA  3.2.  If  *  1«  a  »olutlon  of  (2.1),  then 

g 

*(5>  -  |  +  A,£  1  +  0(cS  as  £  +  0  , 
for  aoae  constant  A,  and  any  t  e  (B.,1).  Here  8.  is  as  In  Lemma  3.1. 


Proof,  (i)  Define 


(  y<5>,  o  <  e  <  i  , 

9(5)  -  { 

^  0  ,  1  <  5  , 


and 


h{£) 


Py<5>, 


0  <  5  <  1  , 


—  /q  r(£,n)Y(n)dn,  1  <  5  , 


so  that  (3.1)  becomes 


9<5)  *  ~  ]n  r(£,n)g(n)dn  +  h(£),  0  <  £  <  "  , 

O  0 


since  the  integrand  is  zero  for  n  >  1,  and  both  sides  are  zero  for  £  >  1.  Setting 

n  »  £u,  we  arrive  at  a  form  suitable  for  the  Mellin  transforms 

g(5)  *  ~  /_  R(u)g(£u)du  +  h(£),  0  <  £  <  •  ,  (3.8a) 

»3 

where 

R(u)  »  r(1,u)  "  -  ~  log  ^  .  (3.8b) 

u 

Because  of  the  plethora  of  notation,  we  emphasize  that  g  is  a  truncated  form  of 
y  “  ^-w/6,  that  h(£)  is  presumably  of  smaller  order  than  g(£)  as  £  ♦  0,  and  that  the 
dependence  of  h  on  r  is  now  implicit. 


1  3. 


(ii)  The  term  h(£)  may  be  estimated  as  follows.  Let  a  e  (— ,— )  be  the  exponent 

2  4 

2. 2)  then  for  £  >  2  (and  hence  u  *  —  when  we  aet  n  •  £u) 

|h<C)|  <  /J|r(5,n)|  |Y(n)|dn 

<  const.  £°  /^^|log  I » 1  | ua  du 


in 


,  .  .a  , 1/F  ,  1  a 

<  const.  £  Jp  log  —  u  du 


<  const.  £  1  log  £  , 
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and  for  1  <  £  <  2  the  integral  is  bounded  by  a  constant.  Since  *  0(£ )  as 


£  +  0,  we  have 


|h(C) |  < 


const.  £, 


0  <  5  <  1  , 


^  const .5  (1+  log  5 ) ,  1  <  5 

It  follows  that  h(s)  exists  and  is  analytic  for  -1  <  a  <  1,  where  s  m  o+it,  and  frost 


(3.6)  that 


h(o+i.)  e  L  ,<-*•,»>  for  all  a  e  (-1,1)  and  all  p*  >  2 
P 


(3.10) 


(ill)  To  find  the  Mellin  transfora  of  the  function  R  in  (3.8),  we  recall  that  our 


original  kernel  in  (2.1)  was 


k(r.,n)  -  1  log 


and  note  from  (121,  p.  192,  that  (with  s  »  o+it) 


Jq  h*  \(1,u)du  «  tan  ~~  ,  -1  <  a  <  1  . 


Setting  u  «  n/C  for  fixed  €  e  (0,«»),  sniltiplylng  both  sides  by  (  ,  and  then  integrat¬ 
ing  with  respect  to  ( ,  we  obtain 

m  -S+1 

/0  h“"  g(C,n)dn  -  tan  §-  .  -1  <  o  <  1  . 

We  Bay  integrate  by  parts  because  n*q( C ,n)  is  0(n°+1  log  1/n)  as  n  +  0,  and  is 


0(n  )  as  n  ♦  •»,  and  there  results 


-  /“  n*qn(C,n)dn  -  tan  ,  -i  <  o  <  i 


Set  n  <*  5u  and  s  «  s-1i  then 


1  im  s-1  ,  ] 1— u  I  .  1  ,  (s-1)» 

-  Jn  u  log  J — Jr*-  du  -  -  tan 


,  0  <  Re  z  <  2  , 


(3.11a> 


that  is. 


Mr)  -  1  tan  i=2ii 


0  <  Re  t  <  2  . 


(3.11b) 


(iv)  We  are  now  in  a  position  to  take  the  Mellin  transform  of  the  integral  equation 
(3.8a).  By  Lemna  2.2,  g(C)  -  0(5®)  as  5*0,  with  a  e  (-j,-^);  therefore  g(s)  exists 

A  A 

and  is  analytic  for  o  >  -a.  We  have  just  shown  that  h(s>  and  R(-s+1)  exist  and  are 
analytic  for  -1  <  o  <  1.  Hence,  if  -a  <  a  <  1,  we  may  apply  the  Mellin  transform  to 


(3.8a)  and  use  the  product  formula  (3.5)  with  A  •  1  and  B  «  0s 

«  1  * 

g(s)  -  —  R(  l-a)g(s)  +  Ms),  -a  <  <J  <  1  , 

/3 


J*  ^ 


•  n-.v.v’av' 


or,  equivalently 


g(s>  -  h Is)  -  (?(s>hfs>. 


where 


Q(s) 


R1 1-s)  _  _ 


-O  <  O  <  1 


tan(s*/2) 


/3 -R  ( 1  -a )  /3(1-s)  +  tan(s*/2) 

The  inversion  formula  (3.4)  qives,  for  all  E  e  (0,»)\{l}, 

g(E)  -  h(  E )  -  “■  e  SQ(s)h(s)ds, 


-a  <  c  <  1 


(3.  12) 


(3.13) 


(3.  14) 


(v)  tet  E  e  (0,1),  so  that  g(E)  »  Y(C).  Since  h(C)  “  0(C)  as  C  ♦  0,  it 

61  i 

suffices  to  show  that  the  right-hand  side  of  (3.14)  behaves  like  A^E  ♦  0(C  )  for  some 
constant  A j  and  any  t  e  (^il).  Me  prove  this  by  moving  the  path  of  integration  as  far 
to  the  left  as  our  knowledge  of  h  allows,  and  taking  the  residue  at  any  pole  of  the 
integrand  between  the  path  in  (3.14)  and  the  new  one.  Hie  details  are  as  follows. 

Given  c  e  (-a, 1)  and  t  e  (8^,1),  define  to  be  the  (positively  directed) 

rectangular  contour  with  corners  at  c  ♦  iM  and  -t  ±  im  for  some  large  M  >  0,  and 
consider  the  contour  integral 


m  A 

27iJr  5  Q(s)h(»)dm  . 


2*1 

n 

The  contribution  of  the  horizontal  parts  <-<  <  cr  <  c,  t  *  ±M)  of  tends  to  zero 

as  M  ♦  «,  because  (a)  |C  *|  -  C”°  with  0  <  C  <  1>  (b>  the  fact  that  tan(s*/2)  *  ±1 
as  t  ♦  t“  with  o  fixed  implies,  in  view  of  (3.13),  that  Q(»)  ■  0(t_1)»  (c)  the 

A 

Riemann-Lebesgue  lemma,  applied  to  the  definition  integral  (3.3)  of  h(s),  shows  that 
h(s)  *0  as  t  *  If  with  a  fixed. 

A 

Since  h(s)  is  analytic  for  -1  <  o  <  1,  we  conclude  from  (3.13)  and  Lemma  3.1  that 
the  only  singularity  of  C  SQ(s)h(s),  with  -1  <  o  <  1,  is  a  simple  pole  at  s  -  -6,. 
Accordingly, 

27T  {/££  -  C"*Q(s)h(s)ds> 

8  .8 

-  C  1  res  Qtsij^g  h(-B1)  -  A^  1 

and  so 


say 


(3.15) 
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(3. 16/ 


P 1  r  ga  _4  4.  * 

g(E>  -  h(E)  -  A^,  ♦  ~  /_,  E  Q(-f4it)h(-i4it)dt  . 

He  bound  this  last  integral  (without  the  factor  E*/2it )  by  means  of  Holder's  inequality, 
noting  that  0(-t+i.)  e  L^^l--,")  because  p(s)  “  0(t_1)  as  t  ♦  i“,  and  that 
h(-i4i.)  e  by  (3.10).  As  we  remarked  earlier,  the  desired  result  now  follow, 

because  g(E)  »  y(E)  for  E  <  1  and  h<€)  *  0(E)  as  E  ♦  0. 

Remarks.  1.  As  was  to  be  expected  from  (3.1)  and  (3.7),  we  cannot  evaluate  the 
coefficient  A1  in  Lemma  3.2 ;  in  (3.15)  we  have  defined  it  in  terms  of  the  global  behavior 
of  y. 

2.  Stokes  conjectured  (see  [1],  pp.  194  and  199)  that  the  profile  of  an  extreme  wave 

is  convex  between  the  crest  and  the  trough;  in  other  words,  that  (i(E)  is  non-increasing 

on  ( 0,“) .  Now  Lemma  3.2  and  analysis  of  the  equation  which  results  from  differentiation 

61  £ 

of  (2.1)  imply  that  E^'(E)  “  A^E  +  )•  Therefore  we  expect  that  A1  <  0,  and  a 

calculation  shows  that  this  corresponds  to  b^AO  in  (1.4).  Numerical  evidence  ([6], 

equation  (2.10)  with  the  value  B  -  0.131  on  p.  776 ;  (14),  Table  1)  suggests  that  b1  >  0 

and  hence  that  A.  <  0. 

28 

3.  The  term  A  E  in  our  next  result  has  a  character  quite  different  from  that  of 

*1 

AfE  ;  instead  of  arisinq  from  (3.7),  it  arises  from  terms  that  are  essentially  sauares  of 
the  first  perturbation,  and  so  A2  li  determined  uniquely  by  A1* 

THBORM  3.3.  If.  ♦  is  a  solution  of  (2.1),  then 

8  28 

♦(E)  -  |  ♦  AfE  1  +  A2E  1  ♦  0(E*)  as  E  ♦  0  , 

for  any  m  C  (28^,8^48).  Here  A^,  8^  and  I  are  as  in  Lemma  3.2,  so  that  i  €  (8^,1), 

and 

^^Xtan  2 

A  ■  •  * ^ ^  0»  1341  A .  e 
2/3(l4B1)  (i'll  1428,) -tan  8^}  ~ 

Proof .  He  shall  use  the  notation  in  the  proof  of  Lensm  3. 2,  and  shall  begin  by 

showing  that  the  previous  estimate  (3.9)  of  h  can  be  sharpened,  by  means  of  the  result  of 

Ismma  3.2,  to 
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0  <  5  <  1 


(3. 17a, 


h(C>  -  +  c2C  +  0(C  ,. 

|h(C)|  <  const.  C_1(1  +  log  i) ,  1  <  ?  ,  (3.17b, 

where  C,  is  a  constant  of  no  interest,  because  ultimately  it  will  drop  cat  of  the 
integral  equation,  and  C2  is  a  constant  that  is  known  in  terms  of  A^ •  With  this 
estimate  in  hand,  we  can  move  the  path  of  integration  in  (3.14)  further  to  the  left  than 
was  possible  before.  Taking  due  account  of  the  poles  of  the  integrand  between  the  path  in 
(3.14)  and  the  new  one,  we  shall  obtain  the  result  of  the  theorem. 

(i)  Let  ?  e  (0,1].  A  straightforward  calculation  gives  the  following. 

<rr)(n)  -  2w(n)  sin{|  +  y(n)) 

2B  6  +i 

■  1  +/3iln)  •  ja^h  +  o(n  )  , 

(FY)(n)  -  /3  y(n)  -  iog{^  J*J  rr)  +  -1-/3  y(n) 

no  ,n 

U  * 

4+86 i+B^  26  6,+i 

- - 1 — £  A,n  +  o(n  )  . 

2(1+Btr 

Since  26 1  and  B^  +  t  are  not  odd  positive  integers.  Lemmas  4.2  and  4.3  yield 

11  -  2fJi  Bi+t 

j /0  r(5,n){(ry)(n)  -  /3  y(fi))dn  -  cQe  +  c2t  +  o(C  )  ,  (3.18) 

where  the  constant  Cg  will  be  a  part  of  C.,  (and  is  therefore  of  no  interest),  and 

(4+86 ,+B*) (tan  B  w)A2 

C2 - 2 -  ' 

6(1+B1)^(1+2B,) 

The  other  tern  in  the  definition  (3.2)  of  P^(n)  is 

2 

4  / 7  rtt.nXPYXn)  -  -  r—  jT  log(1  -  ^)(Fy)(n)dn 


(3.19) 


where  the  integral  is  0(k-^)  as  k  ♦  “,  so  that  the  series  converges  for  5  <  1. 
Equations  (3.18)  and  (3.20)  imply  that 

20  0  +t 

PY<C)  -  ct5  ♦  c25  +  0(5  )  . 

and  thus  prove  (3.17a),  while  (3.17b)  is  a  previous  estimate. 

(ii>  Define 


20. 


h,(5) 


c,C  +  C25 
0  , 


0  <  5  <  1 
1  <  5  . 


and  h2  “  h-hj.  Then,  for  o  >  -1, 


ht(s) 


8+1  8+20, 


while  h2(s)  exists  and  is  analytic  for  -0^  -  l  <  a  <  1,  and  h2(o+i.)  e  L^t  (-**,•)  for 
all  a  e  (-8-0,1)  and  all  p'  >  2.  The  obvious  analytic  (more  precisely,  meromorphic ) 
continuations  of  h  ^  and  h  into  the  strip  of  convergence  of  h2  will  also  be  denoted  by 

A  A 

hj  and  h,  respectively. 

We  can  now  proceed  from  (3.14)  as  in  the  proof  of  Lemau  3.2,  step  (v),  except  that  now 
we  move  the  path  of  integration  to  a  -  -m  tor  any  a  e  (20^,0^+i),  and  collect  the 
residues  of  the  poles  of  5  *Q(a)h(s)  at  s  -  -0  ,  -  1  and  -20^  Since  Q(-1)  *  -1, 


there  results 


6  28 

g(5)  -  h(5)  -  (1,5  1  -  C,5  +  C2Q(-201)5  1 


+  JZm  E"itC(-»+it)h(-si+it)dt  ; 


the  last  integral  is  bounded  as  before  since  h^l-m+i.)  and  h2(-m+i.)  are  both  in 

I*3(-+* ,")•  Substituting  for  h(5)  tram  (3.17a),  we  have 

0  28 

7(5)  -  A,5  ♦  CJ}{O(-201)  +  1)5  +  0(5*)  as  5  +  0  , 

where  C2  is  given  by  (3.19)  and  O("201>  by  (3.13). 
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4. 1 .  Notation 


Let  «0  -  {0,1,2,...  },  II  =  {1,2,3,...  }  and  B  -  {B^.B^Bp. . . .  ),  where  the  numbers 
8^  are  as  In  Lemma  3.1.  We  shall  prove  in  the  Appendix  that  the  B  ^  are  transcendental 
numbers. 

An  elaborate  notation  Is  needed  for  the  exponents  in  various  series  that  we  shall 
construct,  because  (a)  the  non-linearity  in  (2.1),  and  the  presence  of  the  weighting 
function  w  there,  cause  the  exponents  in  the  expansion  of  i|>  to  form  a  fairly 
complicated  sequence,  (b)  further  exponents  enter,  and  finally  depart,  during  the  course  of 
the  construction. 

In  what  follows  (m,,... ,mr+1 )  is  a  multi-index:  mt  e  Nq  for  i  -  1,...,r+1  and 
r  is  any  element  of  M.  Define 

A  -  {.,0,  +...-  mrBr  ♦  mr+12  :  »,  ♦...♦  «r  4  mr+1  >  1}  , 

since  A  is  a  countable  set,  we  may  write  A  -  {^,^'Xj,...  )  with  <  *j+j*  ®*e 

first  few  numbers  are 

Bt,  2B1#  2,  361#  B1  2,  B2,  4B, ,  y  +  2,  0,  ♦  B2,  4,  50^  301  ♦  2,  2B1  +  02 . 

The  following  subsets  of  A  will  be  needed. 


so  that 


and 


A0  -  {.,0,  ♦...♦  mr0r  +  Vl2  :  m,  +  ...♦  -r  >  1>  , 


A  -  A0  u  {2,4,6,...  }  , 


A3  -  {2}  u  {-,0,  ♦...♦  mrBr  ♦  »r+,2  :  m,  +...♦  mr  +  m^,  >  2}  , 


so  that 


A  ■  Aj  u  B  . 

Thus  each  element  of  Ag  involves  at  least  one  of  the  B^ i  we  expect  the  exponents  in  the 
expansion  of  4  to  be  in  A  .  The  exponents  in  A  arise  if,  for  example,  we  combine 

p  3 

three  series  as  follows:  form  the  product  of  two  series  with  exponents  in  A  and  add  a 
series  with  exponents  in  {2,4,6,...  }. 

Truncated  subsets  will  be  denoted  by 


Afl(n)  -  {A  e  A  s  A  <  A  )  -  {X .  e  A  :  j  <  n> 
n  3  B  3  n  3  " 

Ajtn)  ■  {A^  e  Aj  :  3  <  n!  , 

B(n)  -  {8.  e  B  :  8  <  A  )  . 

3  3  n 


4.2.  Bootstrapping  from  y  to 

Ha  now  prove  that  a  given  expansion  of  Y,  with  exponents  in  A  (n),  implies  an 

P 

expansion  of  with  exponents  in  A^(H)  {1,3/5,...  },  where  N  >  n+1.  This  will 

allow  us  to  improve  the  given  approximation  to  y  by  the  method  that  we  used  to  pass  from 

Lews la  3.2  to  Theorem  3.3.  He  shall  prove  in  Theorem  4.5  that,  if  A^(N)  B(N)  »  0,  then 

y  has  an  expansion  with  exponents  in  the  larger  set  A  (N).  Hence,  if  A  B  «  0  (which 

P  3 

is  a  number-theoretic  problem  addressed  in  the  Appendix),  then  y  has  an  infinite 
asymptotic  expansion  with  exponents  in  A^. 

He  begin  by  expanding  the  function  FY  -  /3y  in  the  definition  (3.2)  of  p^  I  then  we 
establish  two  basic  properties  of  the  linear  integral  operator  in  (3.1)  and  in  the  first 
term  of  p^.  The  second  term  of  p^,  has  a  very  simple  expansion.  Combining  these 
results,  we  obtain  the  expansion  of  ,  for  a  given  expansion  of  Y • 

LEMtA  4.1.  Suppose  that,  for  certain  constants  Aj_, 

A 

Y(T,)  ‘  ^A  eAfl(n)  V  +  0(n  1  2*  n  *  °  '  U-1) 

1  P 

where  A  <  u  <  A  and  |i  +  8  e  A  {1,3,5,...  >.  Define  N  “  N(n,u)  by 

A^  -  «a*{Aj  e  A  :  Aj  <  u  +  8^  . 

(Since  A  <  A  +8.  <  ii  +  8.,  we  have  N(n,p)  >  n+1.)  Then 
n+i  n  i  i  " 1  ■ " 

A  p+8 

(FY)(n)  -  /3  Y(n)  -  L  .  B  n  3  +  o(n  )  as  n  ♦  o  ,  (4.2) 

AjCAjlH;  3  — 

where  the  coefficients  B j  depend  only  on  the  A^  _in_  (4.1). 

Proof .  Again  we  abbreviate  (FY)(n)  to  F^(n),  and  similarly  for  other  functions  of 
the  same  kind.  He  have 
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where 


E  In)  -  2w(n)  sin{—  +  y ( n ) }  , 

Y  O 


^  ^  OB 

*<n)  -  (l+n2)  /z(i+bn2)  /z-  i  +  \ 


2k 


k-1 


C^n  >  say,  for  n  <  1 


Define  two  new  operators  H  and  J  by 


H^tn)  -  E^(n)  -1-/3  r(n) 


w(n){cos  Y(n)  -  1  +  /3  [sin  y(n)  -  Y(n)l>  +  (w(n)  -  i}{1  +  /3  Y(n))  ,  (4.3) 


and 


Vn)  {EY  '  °  * 

Let  us  now  restrict  attention  to  those  n,  say  n  <  nQ»  for  which  |j^(n)l  <  “•  Then 


F  (n)  *  IT  t*>  log  (1  +  J  <n>) ] 

Y  dn  Y 

-J?  [n  {JY-IJ?  +  3  Jy  — ■ •  »'  JYS  v*  ' 
'  EY(r,)  -  1  +  d^  {_  2  nJy  +  1  nJy  "*  *  *  }  ' 


whence 


FT<n)  -  /3  Y(n)  -  «Mn)  ♦  (-  j  u2  +  ~  nJ3  -. . .  )  . 


(4.4) 


Although  this  may  seem  a  very  involved  formula,  it  provides  a  good  way  of  expanding 
F^(n)  -  /T  Y(n),  given  an  expansion  of  Y(n). 

Equation  (4.3)  yields 

V”’ " {1  +  Jt  v2k}{‘ ir  + 4?  -••+  /3  [*  si + 


3  y5 

sT  ,} 


+  l  c  n2k{i  +  /I  Y(n)} 
k-1 


Xi  u+ei 

UjeA3(S)  V  +  0(T1  >  ' 


(4.5) 


where  the  constants  Dj  depend  only  on  the  Ai  in  (4.1)  (and  on  the  C^).  Consider  now 
the  remaining  terms  on  the  right  of  (4.4).  Since 

Vn>  mi  Jo {/I T  +  V  * 
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■ 


1 


XH 


fPQ 


m 


equation  (4,5)  ensures  that  the  expansion  of  has  exponents  in  A,  and  hence  that  the 


expansion  of 


(n  +  m  Y  +  H^},  m  >  2  , 


has  exponents  in  A  .  Since  /3  y  +  H  ,  and  hence  J  ,  are  known  with  error  it 

m  m- 1  , .  Y 

follows  that  and  {/3  Y  +  H  }  are  known  with  error  0(n  )  for  m  >  2.  We 


apply  those  two  conclusions  to  the  terns  in  (4.4)  that  follow  H^,  and  the  lemma  is 


proved . 


LEMMA  4.2.  For  K  e  (0,1)  and  p  >  -1, 


*  k(1+p-2k) 


ll  r(e,n)tpdn  -  <  Kp  log  7  +  — 


5  xM+p)2 


k( 1+p-2k) 


if  p  e  {1,3,5,...  }  . 


Proof.  Denote  the  integral  by  and  define 


Cp,  0  <  5  <  1, 


0  ,  1  <  e. 


so  that  ♦  (s)  -  — —  for  o  >  -p 
P  B+p 


I  (€)  -  L  r(C,n)  *  (n)dn 

p  0  p 


(0  <  5  <  i) 


I n  r( 1,u)  ♦  (Su)du 
0  P 


;R(1-s)^p(s)  for  -  min{p,l)  <  <J  <  1  , 


where  we  have  used  the  product  rule,  the  notation  r(1,u)  “  R(u),  and  the  statement  in 


(3.11b)  that  R(1-s)  exists  and  is  analytic  for  -1  <  a  <  1.  By  the  inversion  formula  and 


(3. 11b), 


V5>  “  m  d-  5"8  d8'  “B,in{p'l)  <  c  <  1 


Sciftfi&a  -  •  •  •  -  •  •  •  -  •  •  •  • 


.V.V; V-V 


■-  -V  \  \\  V  *  •  V  v  V  .*  Vv  mN  ,v,N 


V"  '  j,  ’>S  5  * £> 


jeA3(N)n{odd}  Gj5  3  log  X 


*  o{£  >  . 


tan(A  */2) 


Fj  “  3(1+Aj)  Bj*  Gj  “  SXI+X^)  Bj  '  ,4-7) 

but  the  coef f lclentg  Ejj  are  defined  In  terms  of  the  global  behavior  of  Y  (and  not 
merely  in  ter—  of  the  previous  and  Bj ) . 

Remark.  In  (4.6),  the  set  A3(H)  ft  {odd}  is  probably  empty  for  all  N,  but,  if  it 
is  not,  the  exponents  in  it  cause  no  complication  (in  contrast  to  any  exponents  that  there 
may  be  in  the  set  Aj(N)  n  B(N)  to  be  considered  presently),  because  the  logarithmic  terms 
in  (4.6)  will  be  cancelled  when  we  compute  the  corresponding  expansion  of  y. 

Proof.  Recall  from  (3.2)  that 

PY(e>  “  3  ^0  Y(n)>dn  +  ^  J“  r(E,n)(FY)Cn)dn  ,  (4.8) 

and  observe  that,  for  C  e  (0,1), 


2 

/"  r(C,n)<rrKn)d«i  m  -  ~l  1“  logO  -  ^)(nr)(n)dn 


•  2k- 1 

7  l  — j*  (FY)(n)n"2kdn  , 


whsru  th«  last  integral  is  0(k“  >  as  k  ♦  •  Because  (Fy)(n)  ■  0(log  n)  as  n  ♦  •  . 

K 

For  the  first  integral  in  (4.8),  we  apply  Lemma  4.2  to  the  terms  B^n  J  in  (4.2),  and 
Lemma  4.3  to  the  O-term  there,  recalling  that  U  +  B 1  is  not  an  odd  integer.  The  terms  of 
form  c  E  (where  k  6  R  and  cv  ”  0(k”2)  as  k  ♦  •>  that  result  from  the  B^n  ^  in 

*  *  j 

(4.2),  from  the  Oterm  there  and  from  (4.9)  are  all  collected  in  the  first  sum  of  (4.6j 
if  2k -1  <  iwBt  (so  that  k  <  P)  or  in  the  0-term  of  (4.6)  if  2k -1  >  u+fl^  In  the 
second  sum  of  (4.6)  we  have  replaced  A^W)"'.  {odd}  by  AjIHIM  because  (4.7)  shows 
that  Fj  •  0  when  A^  is  an  even  integer. 


4.3.  The  wMBiion  of  Y 


t-— 4.4  provide*  a  generalization  of  step  (i)  in  the  proof  of  Theorem  3.3,  and  we 
are  now  ready  to  generalize  step  (ii>.  To  avoid  a  conceivable  but  highly  improbable 
complication,  we  make 

ASSUMPTION  A(N) .  The  sets  AjIN)  and  B(N)  are  disjoint. 

The  reason  for  this  will  be  explained  after  the  proof  of  Theorem  4.S.  We  have  veri¬ 
fied  Assumption  AIN)  numerically  for  N  <  100,  and  conjecture  that  it  is  true  for  all  N. 
THEORBt  4.5.  If  the  hypotheses  of  Lemma  4.1  and  Assumption  AIN)  hold,  then 

K 

Y<5>  '  ^i.eA.iH)  V  +  0(5  >  52  5  +  0  ' 

3  P 

for  some  v  e  •  Here  a  new  coefficient  Aj,  with  j  euch  that 

X  e  Ag(H)\  Ag(n) ,  is  determined  by  the  previous  coefficients  A^,  with  1  such  that 
Xt  e  it  and  only  if  X^  4  BIN). 

Proof.  Again  we  use  the  notation  in  the  proof  of  Lemma  3.2.  In  view  of  |4.6),  we 
define 


l  V 


2k- 1 


k-1 


+  ^*j®*3 (S>>*  ^  1  +  ^ X 


1^(5) 


j-*  '  ‘•XjeA3IN)n{odd}uj 

if  0  <  C  <  1  , 


OX  Jlog 


0  if  1  <  C, 


and  ha  “  h-hj.  Accordingly,  hjl?)  -  P^(C)  -  h1!5)  it  C  <  1,  and  h2(C)  - 

0I5-1  log  K)  as  C  ♦  ",  by  (3.9).  The  smallest  exponent  in  h^lO  is  1>  hence,  for 

c  >  -1, 

h1<S>  *  a+2k-1  *  ^X3eA3IN)\*  s+X^  +  ^eAjINJfl  {odd}  (g+x  }2  ' 

A  A 

while  hjls)  exists  and  is  analytic  for  -W-6^  <  o  <  1,  and  hjlo+i*)  e  L^,  (-*,")  for 
all  o  e  l-u-B^.I)  and  all  p'  >  2.  The  obvious  analytic  (more  precisely,  meromorphic) 

A  A  A 

continuations  of  h^  and  h  into  the  strip  of  convergence  of  hj  will  also  be  denoted  by 

A  A 

hj  and  h,  respectively. 
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Let  £  e  (0,1).  we  proceed  from  (3.14)  and  move  the  path  of  integration  to  a 


-v 


for  any  v  e  (X  ,u+B.)»  the  remarks,  in  the  proof  of  Lemma  3.2,  regarding  the  contribution 
N  i 

— g  * 

of  paths  t  -  +H  remain  valid  when  the  t  there  is  replaced  by  v,  and  £  Q(s)h(s)  has 

no  singularities  for  a  e  (-u-6.,-A  ).  There  are  four  sets  of  poles  between  the  original 

I  N 

path  of  integration  (a  -  c)  and  the  new  one  (a  «  -v)»  we  consider  one  set  at  a  time. 

(a)  The  points  s  -  1-2k,  k  B  {1,2,  ..,!>},  are  simple  poles  of  h^  and,  since 
Q(1-2k)  -  -1, 

}  residues  «  -  £  ££2k  1 

(i>  k-1 

(b)  The  points  s  -  -A^ ,  where  A^  e  A^tlDS.  ■,  are  simple  poles  of  h1  and  are  not 

poles  of  Q  by  Assumption  A(H),  so  that 

I  residues  -  L  (M)V-  FjQ(-Xj)£  3  . 

(b)  j  3 

(c)  The  points  s  -  -A^,  where  X^  e  A3(H)  n  (odd),  are  double  poles  of  ht,  and 

l  residue.  -  -  lX  eK  ,v)n{oAi)  6jU  3  109  £  +  (1+Xj)5  3}  * 

(c)  j  3 

(d)  The  points  s  -  -A^,  where  X^  B  B(H),  are  simple  poles  of  9  and  distinct  from 

those  of  h  by  Assumption  A(H),  so  that 

1  X 

l  residues  •  L  H..E  J  • 

(d)  Aj8B,1,)  3 

where 


-  res  OC * > | »— X  *  (4.10) 

Since  £  B  (0,1),  we  have  g(£>  -y<£)  and  h(£)  -  > »  equation  (3.14)  yields 

f(£)  -  PT(£)  +  l  residues  +  /_vti»  5  “c(s)h(s)ds  . 

2k- 1 

We  substitute  for  P  (£)  from  (4.6)  and  for  the  residues,  noting  that  the  terms  E.  £  , 

Y  * 

and  those  containing  log(1/£),  in  (4.6)  are  cancelled  by  the  residues;  there  results 


*<*>  “  ^XjeA3(H)X«  V14Q("VK  *  ’  ^j8A3(H)n{odd}  V,+V? 

X.  Fv  -  V*. 

+  VeB(H)  V  1  +  5  Q(-v+it)h(-v+it)dt  "  0<e  »  ' 


(4.11) 


where  Q(-v+i.)  8  L  (-“,«•),  because  Q(-v+it)  »  0(t  )  as  t  ♦  +»,  and 

3/2 

*  V 

h(-v+i.)  e  ,  so  that  the  integral  (without  the  factor  5  /2w)  is  bounded,  by 

Holder 1 s  inequality • 

He  claim  that  all  the  exponents  1.  in  (4.11)  belong  to  A  (N).  Indeed,  if 

3  P 

X.  e  A\A„,  then  1.  e  (2,4,6,...  },  and  this  is  impossible  in  (4.11)  because  there 
3  P  3 

\  4  *  or  X  e  {odd}  or  X  e  B.  Since  in  (4.11)  each  X  <  X  ,  it  follows  that 

3  3  j  3  ■ 

X.  e  A  (N). 

X3 

Comparison  with  (4.1)  shows  that  in  (4.11)  the  coefficient  of  C  must  equal  A., 

X, 

whenever  X^  e  A^(n).  Consider  the  coefficient  o(  5  J  when  X^  8  Ag (N)N.  Ag (n) j  if 
X^  ^  B(N),  then  (4.7)  and  Lmna  4.1  show  that  the  coefficient  is  determined  by  the 
previous  A^j  if  X^  e  B(N) ,  then  (4.10)  shows  that  the  coefficient  is  not  so  determined. 
Of  course,  we  relabel  this  new  coefficient  Aj  in  either  case. 

Remark.  Assume  that  the  hypotheses  of  leau  4.1  hold  but  that  Assumption  A(H)  is 
false,  for  some  particular  n  and  N(n,li).  Then  there  exists  an  exponent 

A 

X^  e  A^(N)  n  B(N) ,  and  Q(a)h(s)  has  a  double  pole  at  s  »  -X^.  mis  causes  the  term 

with  j  -  m  in  (4. 11)  to  be  replaced  by 

X 

TJ.  *(a+b  log  £)  , 

where  a,  b  are  constants  and  b  j*  0.  At  the  Kth  stage  (and  perhaps  earlier),  the  hypo¬ 
thesis  (4.1),  with  n  replaced  by  M(n,M),  must  be  modified  to  Include  the  logarithmic 
termi  at  still  later  stages,  higher  and  higher  powers  of  logarithms  accrue  in  the  process 
of  expanding  FT  -  /5y. 

COROLLARY  4.6.  _If  ♦  is  a  solution  of  (2.1)  and  Aj  n  B  *  I,  then 


♦<s>  -  i  * )»,«,  V  1  -  '  * 0 

3  P 


(4.12) 


in  the  sense  of  an  asymptotic  expansion  > 


for  all  new. 


♦<*>  "  ?  ‘  Ix^fn)  V  "  0<?  B>  22  5  +  0  ' 


(4.13) 


COROUARY  4.7.  _If  n  B  ■  flf,  then  the  complex  co-ordinate  z  -  x+i y  is  rela 
the  ccwplex  potential  x  “  •♦if  by 

z(X>  ~  -ig"1/3(|  iX)2/3(l  +  );x  @A  ^(ix)  3)  as  x  ♦  0  , 

j  6 

where  the  constant  g  Is  the  gravitational  acceleration,  arg(ix)  e  f-*/2,x/2) ,  and  the 

coefficients  Kj  are  real.  If  the  set  {l} u  B  is  linearly  Independent  over  the 

rationale,  then  the  sum  in  (4.13)  la  only  over  exponents  in 

A-  „«  {n.0,  +...♦  m  6  :  m.  +. ..+  m  >1}  t 

6,0  11  r  r  1  r 

in  other  words,  KU  •  0  whenever  A.  e  A.\A. 

1  J  — — .  j  D  Pf0 

Proof.  On  the  inage  f  «  0  of  the  free  boundary,  we  have 

«• 

£  -  -a4{1  ♦  l  c  *2k),  a  >  0  ,  (4.14) 

k-1 

where  a  and  are  constants,  and  the  series  converges  for  sufficiently  snail  values  of 

|4|.  Define  6(4)  -  f(£(4))s  it  is  a  basic  hypothesis  in  the  derivation  of  (1.1)  that  x 
and  6  are  odd  functions  of  •,  while  y  is  even.  Accordingly,  (4.12)  and  (4.14)  isiply 


that 


•(♦> - sgn  4{^  ♦  lx  gA  B^|*|  J>  as  ♦  ♦  0 


(4.15) 


for  certain  coefficients  Bj  (which  are  not  those  in  (4.2)).  Moreover,  it  follows  froet 
the  boundary  condition  of  the  basic  problen  (or,  equivalently,  from  (1.1),  frost  the  fact 
that  6(4)  is  the  boundary  value  of  Is  log(dz/dx),  and  frost  suitable  choice  of  the 
additive  constant  in  Re  log(ds/dx>)  that 

If  -  (-3g  jj  sin  6)"1/3ei8,#>,  *  e  [-c,c]\{0),  f  -  0  ,  (4.16) 

for  sosm  constant  c  >  0.  using  (4.15)  in  this  formula,  we  obtain  (4.13)  for  z(4)  (that 

is,  for  arg(ix)  “  -»/2  or  »/2 ) i  that  the  are  real  follows  from  the  symmetry . 

To  prove  that  the  asyaptotic  series  for  z(x>  in  (4.13)  is  the  only  appropriate 
extension  (into  the  half-plane  T  <  0)  of  the  series  for  z(4),  we  let  h  denote  the 
difference  between  two  such  extensions  to  finitely  many  terms.  It  then  suffices  to  prove 
the  following!  if  h(x)  is  analytic  and  bounded  in  D  «  (x  s  0  <  |x |  <  c, 

-  ^  <  argdx)  4  ^},  and  if,  for  some  constant  V  >  0,  we  have  h(4)  «  o(|4|U)  as  4  ♦  0 

(with  arg(i6)  “  -*/2  or  w/2),  then  h(x)  ”  0( | X I ^ >  *■  X  ♦  0  in  this  follows 
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from  application  of  tha  Fhragmfcn-LindelSf  theorem  ((11],  p.  176)  to  the  function  f 

defined  by  f(x)  “  *  euitable  auxiliary  function  is  «(x>  “  exp{-(ix)  . 

It  remains  to  prove  that,  if  the  eat  { 1}  u  B  is  linearly  independent  over  the 

rationale,  then  *  0  whenever  _  e  A^\  Ag  Suppose  that 

a  a  . 

s(«)  -  -ig"V3(|  i*)2/3{1  ♦  b,(i*>  1  ♦...+  b  (i*)  n  ♦  k(i*)  ♦  0{  (•  |  >W)  , 

2  i  n 

where  0  <  a.  <...<  a  <  X  <  u,  each  a  e  A  and  1  6  A0\  A-  _.  We  know,  fro*  our 
I  n  ]  PfQ  PDpU 

construction  by  way  of  (4.16),  that  this  approximation  may  be  differentiated  term  by  term, 
and  that  the  exact  function  *(•)  satiafies  the  boundary  condition  of  the  basic  problem. 
Then  a  slight  variant  of  a  calculation  by  Norman  ([9],  p.  262)  shows  that,  under  the 
foregoing  hypothesis,  the  coefficient  k  -  0.  Repeated  application  of  this  argument  provea 
the  result. 
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Appendix .  On  the  numbers  6  ^ . 

We  recall  that  0^  denotes  the  j-th  positive  root  of  /3(1+B)  ”  tan(B*/2)»  Its 
simplest  properties  were  noted  in  Lenna  3.1.  We  now  probe  a  little  deeper. 

LEMMA  A. 1 .  Each  number  8^  is  transcendental. 

Proof.  We  abbreviate  8^  to  8  for  any  fixed  j  e  N.  It  was  shown  by  Grant  ([4], 
p.  260)/  albeit  somewhat  tersely,  that  8  is  irrational.  To  prove  it  transcendental, 
assume  the  contrary:  that  8  is  an  algebraic  irrational  number.  The  equation  of  which 


0  is  a  root  may  be  written 


1*8  1+i/3(1+6) 

e  -  - 

1-i/3{1+6) 


Here  the  right-hand  side  is  an  algebraic  number  (because  8  is  one,  and  the  algebraic 

B 

numbers  form  a  field),  so,  therefore,  is  the  left-hand  side.  In  other  wordB,  (-1)  is  an 
algebraic  number.  But  this  contradicts  the  Gelfond-Schneider  theorem  ([10],  p.  76),  which 
states  that  ab  is  transcendental  whenever  a  is  algebraic  (and  neither  0  nor  1) 
and  b  is  algebraic  and  irrational. 

Remark .  In  Corollary  4.6  we  assumed  that  the  sets 

A3  -  {2}  0  (m^  ♦...+  *r8r  +  *^,2  ,  m.j  »r  +  Vi  >  2} 


B  ■  } 

are  disjoint.  It  is  clear  that  the  truth  of  this  assumption  would  be  implied  by  the  truth 
of  the  following,  very  natural 

Conjecture .  The  set  {l}  u  B  is  linearly  Independent  over  the  rationale. 

(That  is,  if  x1,...,x)(  are  distinct  elements  of  {1}  u  B,  and  a  are 

rational  numbers,  then  \  .ax  -  0  implies  that  a  -...«  a  «  0.)  Unfortunately,  our 

n*i  n  n  i  k 

only  result  in  this  direction  is  the  following. 

THEOREM  A.2.  _lf_  j  ¥  k,  then  the  set  {I.B^.B^}  is  linearly  independent  over  the 
rationale. 

Proof .  Assume  the  contrary:  then  there  exist  integers  p,  q  and  r,  not  all  zero, 
such  that  p0  +  qfl  -  r.  Neither  p  nor  q  can  be  zero;  we  may  suppose  that  r  >  0: 

J  * 


then  at  least  one  of  p  and  q  must  be  positive,  say  p  >  0.  Mow 


,  «vr  i*r 
(-1)  -  e 


lir 0.  p  1*0,  q  l+i/Id+BJ  p  1+i/3d+BJ  q 

(•  ’)  (e  k)  ■  1— - M  i— - H 

l-i/Td+8  l  i-i/T(  1+6.  ) 
j  * 


by  (A. 1)»  since  »  <r-pB^)/q,  we  have 

r-pS  .  q  _  r-pB.  q 

(-1  )r  {l-i/3  ( 1+0  , )  }**{  1-i/3  ( 1  +  - *>)  -  (1+i/3(1+B.)}p{l+i/3<1  + - *>}  .  (A. 2) 

3  q  J  H 

Suppose  that  q  >  Oj  expanding  both  sides  of  (A. 2),  and  re-arranging  the  result,  we  obtain 

i/jp,^)*^)  , 

where  P1  and  P2  are  polynomials  with  (real)  rational  coefficients,  and  0^  is  real. 
Hence  each  polynomial  is  zero,  which  makes  8^  an  algebraic  number  and  thus  contradicts 
lesmw  A.  1. 


If  q  <  0,  we  re-write  (A. 2)  to  have  positive  exponents  p  and  -q  on  both  sides, 
and  argue  as  before. 
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